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ABSTRACT 

Lopsidedness of the gaseous disk of spiral galaxies is a common phenomenon in 
disk morphology, profile and kinematics. Simultaneously, the asymmetry of a galaxy's 
stellar disk, in combination with other morphological parameters, has seen extensive 
use as an indication of recent merger or interaction in galaxy samples. Quantified 
morphology of stellar spiral disks is one avenue to determine the merger rate over much 
of the age of the Universe. In this paper, we measure the quantitative morphology 
parameters for the H i column density maps from the Westerbork observations of 
neutral Hydrogen in Irregular and SPiral galaxies (WHISP). These are Concentration, 
Asymmetry, Smoothness, Gini, M20, and one addition of our own, the Gini parameter 
of the second order moment {Gm)- Our aim is to determine if lopsided or interacting 
disks can be identified with these parameters. Our sample of 141 Hi maps have all 
previous classifications on their lopsidedness and interaction. 

We find that the Asymmetry, M20 and our new Gm parameter correlate only 
weakly with the previous morphological lopsidedness quantification. These three pa- 
rameters may be used to compute a probability that an Hi disk is morphologically 
lopsided but not unequivocally to determine it. However, we do find that that the 
question whether or not an H i disk is interacting can be settled well using morpholog- 
ical parameters. Parameter cuts from the literature do not translate from ultraviolet 
to H I directly but new selection criteria using combinations of Asymmetry and M20 
or Concentration and M20, work very well. 

We suggest that future all-sky H i surveys may use these parameters of the column 
density maps to determine the merger fraction and hence rate in the local Universe 
with a high degree of accuracy. 

Key words: galaxies: fundamental parameters galaxies: spiral galaxies: structure 
galaxies: interactions galaxies: kinematics and dynamics 



1 INTRODUCTION 

In the study of the 21 cm emission of atomic hydrogen (H l) 
from nearby galaxies it was noted early on that many ap- 
pear to be not symmetric. This phe nomenon was terme d 
"lopsidedness" of t he Hi morphology (|Baldwin et al.|[l980l ). 
iRichter fc Sancisil l|l994l ) find that half of galaxy disks are 
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lopsided. A similar deviation from the axi-symm etry in 
the st ellar disks of galaxies was noted by Rix & Zaritskvl 
(| 19951 ) and IZaritskv fc Rb3 (| 19971 ). The Hi line profile of 
half the population of galaxies, also shows a clear devia- 
tion from symmetry on either side of the s ys temic velocity 
( Havnes et all 1 199^ : [Matthews et al.l Il998l ). I Swaters et al.l 
1 1999h report a third kind of lopsidedness, a deviation 
from axi-symmetry in the position-velocity diagram, termed 
kinematic lopsidedness. The fraction of lopsided galaxies 
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may depend on environment; e.g., the Eridanus Galaxy 
^roup counts twice as many lopsided galaxies as the field 



group counts twice as many lopsided galaxies as the neld 
ijAngiras et"alll2006h. but the U rsa Major group is similar 
to the field ( Angiras et al.l [20071 ). Lopsided ness seems to be 
the strongest in the outer regions of a disk (|joelll999l ). 

Disk lo psid edness may be the prod uct of tidal in- 
terac tions (|jod [1997), minor mergers (jZaritskv fc Rixl 
[1993), as ymmetric accretion of fresh gas from the cos- 
mic web (iBournaud et aL 2005al ). intergalactic gas ram- 
pressure (iMapelli et al.ll2009l ~ or an offset between the disk 
and d ark matter halo ("disk s loshing", iLevine fc Spark3 
1 19981 : iNoordermeer et al.l l200lh . Alternatively, the phe- 
nomenon may be at tributed equally to most of these causes 
ijMapelh et al.ll2009l ). In the case of a stellar disk, lopsided- 
ness may h ave an internal cause , such as a dynam ical insta- 
bility fe.g.. lLovelace et al.|[T999l : lDurv et al.ll2008l '). 

Initially, the main way to identify lopsidedness has been 
a visual inspection of the Hi colu mn density map, line pro- 
file or velocity field of a galaxv dRichter fc Sancis 1 1 19941: 
Havnes et al] Il998l: [Matthews et al.l Il998l : ISwaters et all 



20021 : INoordermeer et al.ll2005bl ). In addition, lopsidedness 



can be quanti fied using a Fourier decomposition of either th e 



stellar image (jZaritskv fc Rixl 



1997l:lBournaud et al.ll2005al ). 



the Hi colum n density map (Angiras et aL[2006l. 20071) . or 



veloc ity field l|Schoeninakers et al. ll997l : lTrachternach et all 
I2OO8I ). However, the Fourier analysis has only been per- 
formed on small samples of H I observations or larger sam- 
ples of op tical ones. We r efer t he reader to the in-depth 
review by |jog fc Combed (|2009l ) on the lopsidedness phe- 
nomenon. 

In parallel with the line of investigation into galaxy 
lopsidedness, a considerable observational effort has gone 
into morphological tracers of interaction over cosmologi- 
cal times. These studies use certain quantifiable morpho- 
logical parameters of restframe-ultraviolet images of distant 
and nearby galaxies to e s timate the merger ra t e of galax- 



T ge 

ies j Abraham et all 1 19941 : IConselice et al.|[200ol : iLotz et al] 



l2004 r Two sets of parameterisations ha ve emerged, the 
Concentration- Asymme try-Smoothness by IConselicg (|2003l ) 
and the Gini-M2o by iLotz et al] (|2004l ). These parame- 
terisations of galaxy morphology have now been applied 
on every deep multi-wavelength Hubble field to deter- 
merger rates ; in the Hub ble Ultra Deep Fiel d by 
Pirzkal ot al. (200^, GOO DS bv iBundv eraP ('2005') and 



Oj). 

Rayindranath ct al] (|2006|) . COS MOS by [Scarlata et al 



and IConselice et all ((iooi), GE MS bv IJogee et al.l 
and the extended G roth strip bv lLotz et al.l ( 2008al ) 



Jogee et al.l 



and Conselice et al.l (|2008|). a s w ell as l ocal reference sam- 



ples ([Conselice 



20031: 
it aLl ls 



unoz-Mateos et 
of galaxy appearance are sensitive 



.Lotz ct al 
20091 ) 



2004; Hond o et all l2007l : 

The different parameterisations 
to differen t stages of 
an i nteraction and different interaction types (|Lotz et al.l 
l2008 b: ConscUcG 2009; Lotz et al. 2010a, b) but are very suc- 
cessful in estimating the galaxy merger fraction and rate 
over much of the age of the Universe. 

Yet, to date, these studies have been constrained 
mostly to restframe ultraviolet and optical because these are 
the wavelengths where interaction-induced star-formation 
produces high-surface brightness features in galaxies with 
clearly disturbed morphology at wavelengths where the 

Hubble Space Telescope can reasonably obs erve them. 

In the previous papers in this series (|Holwerda et al.l 



I2OO9I , l2011b| [3). we have shown that a description of the 
Hi morphology using these parameters is as sensitive, if 
not better than, any of the star-formation dominated wave- 
lengths to the effects of intera ctions. Hence, the future 
Square Kilometer Array (SKA; ICariUi fc Rawlingd |2004| ) 
and its precursor radio tele scopes. South Afric a 's Karoo Ar- 
ray Telescope (Mee rKAT; [Booth et al] [20091 : [Jonasi [2007[ : 
[de Blok et"al] [2009[) . and t he Australian SKA Pathfind er 
(ASK AP; [Johnstorj [2007[ : [Johnston et al. [ [20071 . [2008a[ [bl. 
[2009t ). provide an opportunity to explore the lopsidedness 
phenomenon as well as interactions using 21 cm line emis- 
sion (Hi) of thousands of galaxies. However, to do so, auto- 
mated parameterisations of the H I maps are needed and the 
relation between the lopsidedness phenomena and the above 
parameter space will need to be explored. If, for instance a 
parameter space can be identified for interacting galaxies, 
a merger fraction for a given cosmic volume can be mea- 
sured. Combined with an estim ate of the timescale a m erger 
spends in this parameter space ([Holwerda et al. 2011d), on e 
can then estimate the merger rate ( Holwerda et al. 201 la[ ). 
Our ultimate goal is to simplify the selection of subsamples 
in the upcoming large H I surveys using existing morpholog- 
ical parameters. 

In this paper, we compare the CAS and Gini-M2o pa- 
rameters as determined in the Hi column density maps 
of 141 galaxies, to the lopsided ness qualifica ti on an d in- 
teraction determinations from iSwaters et al.l (|2002l ) and 
[Noordermeer et al. [ ([2005b[ ). In 33 we define the two main 
concepts, in ^Swe briefly discuss the morphological param- 
eters and present a new additional parameter. In <21 we de- 
scribe the radio data used. Our results are presented in ^ 
with our conclusions in ij6] and a brief outlook in ^Zl 



2 LOPSIDEDNESS AND ASYMMETRY 

In the following discussion, we find it useful to define the 
terms lopsidedness and asymmetry as they are often used 
interchangeably and the difference is subtle. 
Lopsidedness is a comparison of OEi-symmetry, the level of 
symmetry of a galaxy image, line profile o r velocity field 
when mirrored over an axi s (minor or major) ([Baldwin et al.[ 
[l980l : [Swaters et"al][ 19991 ). Quantified definitions of lopsid- 
edness are when a disk displays an m = 1 global spatial 
offset (m is the azimuthal wavenumber in a spatial Fourier 
decomposition) or the cos (</!>) distribution {(p is the azimuthal 
angle) is non-axisymmetric. In the case of the WHISP sam- 
ple, there are qualitative est i mates of a disk's lopsidedness 
available from [Swaters et "al] ([20021 ) and [Noordermeer et al.[ 
(,2005bl ). 

Asymmet r y, as defined by [Abraham et al.l ([19941 ) and 
[Conselicg ( 2003[ ) . is the point-symmetry of an object; the 
level of symmetry when the object is rotated 180° around 
its centre. 



3 MORPHOLOGICAL PARAMETERS 

The morphological parameters we compute over the H I 
column density m aps are C o ncentr ation, Asymmetry and 
Smoothness from [Conselicd ([2003[ ). M20 and Gini from 
[Lotz et al] ([2004[ ) . and a single addition of our own; the Gini 
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parameter of the second order moment of the hght, Gm- 
W e describe our implemen tation of the existing parameters 
in lHolwerda et al] (|2011bl l3) and below. The relevant input 
parameters are the central position of the galaxy (xc, j/c), 
and a definition of the area over which these parameters 
are computed. The Gini parameter only requires the defini- 
tion of the area and not the central position, making it less 
sensitive to input error. We obtained uncertainty estimates 
from a Monte-Carlo run, varying the central position of each 
galaxy and a separate run randomly redistributing the pixel 
values in the galaxy-area. 



3.1 CAS 

CAS refers to the now com monly used Co ncentration- 
Asymmetry-Smoothness space (jConselicell2003l ) for morpho- 
logical analysis of distant galaxies. Concentration of the 
light, symmetry around the centre and smoothness as an 

indication of substructure. 

Concentration is defined bv lBershadv et all l|200d ) as: 



C = 5 log(r8o/7-2o) 



(1) 



with r/ as the radius containing pe rcentage / of the light of 
the galaxy (see definitions of r/ in lBertin fc Arnoutslll99a : 
lHolwerda|[2005l ). 

The asymmetry is defined as the level of pomt-, (or 
rotational-) s ymmetr y around the c entre of the galaxy 
l| Abraham et aL, 1994 : iConselicdl2003l '): 



S,,i|/(i,j)| 



(2) 



where is the value of the pixel at the position i,j in 

the image, and /i8o(*,i) is the pixel at position [i,j] in the 
galaxy's image, after it was rotated 180° around the centre 

of the galaxy. 

Inspired by the "unsharp mas king" techn i que (MalinI 

19781). Smoothn ess is defined by iTakamival l|l99£ ) and 



Conselicd (|2003l ) as: 

^^,j\Iii,j)-Is{^,j)\ 



S = 



-id\Hi,j)\ 



(3) 



where Is{i,j) is the same pixel in a smoothed image. What 
type of smoothing is used has changed over the years. We 
chose a fixed 5" Gaussian smoothing kernel for simplicity. 



3.2 Gini and M20 

[Abraham et~aLl (|2003l ) and lLotz etal] l|2004 ) introduce the 
Gini parameter to quantify the distribution of flux over the 
pixels in an image. They use the following definition: 



G = 



In{n - 1) 



E,(2i-n-l)/,, 



(4) 



li is the value of pixel i in an ordered list of the pixels, n 
is the number of pixels in the image, and I is the mean 
pixel value in the image. We chose this definition as it is the 
computationally least expensive. The Gini parameter is an 
indicatio n of equality in a distributio n (initially an economic 
indicator [Ginil(l912l : lYitzhakilll99lf l. with G=0 the perfect 
equality (all pixels have the same intensity) and G=l perfect 
inequality (all the intensity is in a single pixel) . Its behaviour 



is therefore in between that of a structural measure and 
con centration. 

iLotz et al] (|2004l ) also introduced a new way to param- 
eterize the extent of the light in a galaxy image. They define 
the spatial second order moment as the product of the inten- 
sity with the square of the projected distance to the centre 
of the galaxy. This gives more weight to emission further 
out in the disk. It is sensitive to substructures such as spiral 
arms and star-forming regions but insensitive if these are 
distributed symmetrically or not. 

The second order moment of a pixel i is defined as: 



M, = hx [(x - x^Y + [y ^ y^Y], 



(5) 



where [x, y] is the position of a pixel with intensity value h 
in the image and [2;c,2/c] is the central pixel position of the 
galaxy in the Hi surface density map. 

The total second order moment of the image is given 

by: 

Mtot = T.iMi = T.Ii[{x, - xcf + {y^ - ycf]. (6) 

iLotz et al.l (|2004 ) use the relative contribution of the 
brightest 20% of the pixels to the second order moment as 
a measure of disturbance of a galaxy: 



M20 = 



V Mtot J 



for E,7, < Q.2Itot 



(7) 



The M20 parameter is sensitive to bright regions in the out- 
skirts of disks and thus higher values can be expected in 
galaxy images (in the optical and UV) with star-forming 
outer regions as well as those images of strongly interacting 
disks. 



3.3 Gini of the second order moment (Gm) 

Instead of using the intensity of pixels, we can define a Gini 
parameter for the second order moment of each pixel by 
substituting Mi (equation [5]) for U in equation U) 



Gm ~ 



Mn(n - 



-E,(2i 



1)M., 



(8) 



This is our contribution to the parameter space to pro- 
vide an additional handle to characterise lopsidedness and 
interaction level. Our reasoning was that the Gini param- 
eter has the added benefit of using the combined shape of 
the fiux distribution curve (all th e information in the im - 
age), rather than just a fraction. In lHolwerda et all (|201ld ). 
we found hints that M20 may be not sensitive enough to in- 
teraction signature while Asymmetry is sensitive to other ef- 
fects as well. A similar conclusion was reached bv lLotz et al.l 
(|2008d ). hence Gm is an attempt to define a single param- 
eter to detect interaction using all the information on the 
second order moment. 



4 WHISP DATA 

The data we use are the Hi column density maps from 
the Westerbork Hi Survey of Irreg u lar an d SPiral galax- 



ies ( WHISP. 



van der Hulst et al.l . _ _ . . . . .- 

20021: [Swaters et al.l |2002|: ISwaters fc Balcellj 1200^ 



2001 : van der HulstI 



Noordermeer et "al] 2005bl ). WHISP is a survey of the 
neutral hydrogen component in spiral and irregular galaxies 
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with the Westerbork Synthesis Radio Telescope (WSRT). 
It has mapped the distribution and velocity structure 
of Hi in several hundreds of nearby galaxies, increasing 
the number of Hi observations of galaxies by an order 
of magnitude. The WHISP project provides a uniform 
database of datacubes, zeroth-order and velocity maps. Its 
focus has been on the structure of the dark matter halo as 
a function of Hubble type, the TuUy-Fisher relation and the 
dark matter content of dwarf galaxies. 

The WHISP observation targets wer e selected fro m the 
Uppsala General Catalogue of Galaxies (|Nilsonll 19731 ). with 
blue major diameters > 2'.0, declination (B1950) S > 20° 
and flux densities at 21-cm larger than 100 mjy, later low- 
ered to 20 mJy Observation times were typically 12 hours of 
integration. The galaxies satisfying these selection criteria 
generally have redshifts less than 2 0000 km/s (z < 0.07). A 
further prerequisite was th at either [Swaters et al. I (|2002) or 
iNoordermeer et al.l (|2005bl ) classified both the level of the 
galaxy's lopsidedness and whether or not it is interacting. 

The WHISP data was retrieved from the 
"Westerbork on the web" project at ASTRON 
(http://www.astron.nl/wow/). We use the column 
density maps with the highest resolution available (~12" 
X 12"/sin(5)). The positions and b asic H I information 
(mas ses and diameters e t c.) are from ISwaters et al] (|2002l ) 
and INoordermeer et al.l l|2005al lbh. We used the central 
position (xc, yc) as input for the parameters and the radius 
of the Hi disk (Rhi) to cut out a stamp of the disk before 
computation (a stamp was set at 7 x Rhi)- The computed 
morphological parameters are in Tables Al and A2 in 
Appendix A {electronic edition only). 



5 RESULTS 



The samples fr om ISwaters et al.l (|2002l ) and 
INoordermeer et al.l (|2005bl ) both have visual classifica- 
tions of a galaxy's lopsidedness and whether or not it is 
interacting. Lopsidedness was determined for each galaxy's 
morphology in the Hi column density map, its profile and 
velocity map. These are the morphological, profile and 
kinematic lopsidedness respectively. These classification 
were done by a single observer by visual inspection and 
he nce carry some risk of observer bias. The classifications 
by ISwaters et all (j200^) for lopsidedness are: not-, weak-, 
an d strong lopsidedne s s. The lopsidedness classification 
by INoordermeer et al.l l|2005bl ) is a little more nuanced 
with no, mildy-, moderately-, and severely lopsided. Both 
authors classify morphological, profile as well as kinematic 
lopsidedness with their respective qualifiers. To unify the 
two cl assification schemes, we re-assigned the lSwaters et al] 
||2002|) classification s to lopsidedness categories of the 
INoorderm eer ct a l] (|2005bl ) classification: weak is equiva- 
lent to mild and strong to severely. In the following section 
we use the Noordermee r et al. l opsidedness scale. 

INoordermeer et al] l|2005bl ) also gives an estirnate on 
whether the galaxy is interacting and lSwaters et al] (|2002| ) 
lists the five galaxies in their sample of 74 dwarfs that are 
in an active interaction. This fraction (5/74) might be an 
under-estimate and may not account for galaxies that are 
only mildly interacting. 

In the following section (15. If) , we compare the visual 



classifications of lopsidedness to our morphological parame- 
ters to determine if lopsidedness can be quantified with our 
morphological parameters. In the next section (|5.2p . we ex- 
plore the distribution of interacting galaxies in our parame- 
ter space. These visual classifications are subject to possible 
observer bias but the aim here is to identify the parts of 
the morphological parameter space described in section [3] 
that hold the majority of lopsided or interacting galaxies. 
Appendix A (electromc edition only) lists the morphological 
parameters of all the galaxies in the Swaters et al. (2002) 
and Noordermeer et al. (2005) sample in Tables Al and A2 
respectively. 

5.1 Lopsidedness 

Figure [T] shows the histograms of the above six parameters 
(C, A, S, G, M20, and Gm) for lopsided (any strength) 
and non-lopsided g alax y morphology acco r ding t o either 
ISwaters et al] l|2002i ) or INoordermeer et all (|2005bl ). Based 
on the definition of lopsidedness, we expected the Asymme- 
try, M20, and possibly the Gm parameters to show a differ- 
ence between the two populations. We observe a difference 
in the median in these parameters with the lopsided galax- 
ies showing higher values for Asymmetry and Gm and lower 
values for M20 (Figure [T} . The Gini parameter shows a shift 
in the median value for the lopsided galaxies as well. 

Figure [2] shows the parameter space of C, A, G, M20, 
and Gm with the different lopsidedness classifications (no, 
mild, moderate and seve rely). As men ti oned above, weak 
and strong according to .Swaters et al.l (|2002l ) are plotted 
as mild and severely respectively. There is no clear part of 
parameter space where one could identify only, for instance, 
the severely (strongly) lopsided galaxies. 

While the distributions of morphological values are dif- 
ferent between lopsided and not lopsided galaxies, there is 
no clear cut in morphological parameters to discern be- 
tween the two or separate out weakly and strongly lop- 
sided galaxies. As it stands, the distributions in Figure 
[T] could be used to compute a probability that a galaxy 
is lopsided, but this would still ha ve to be followed up 
with a visual inspection like th ose in ISwaters et al.l (|2002l ) 
and INoordermeer et al.' ('2005b') or a Fourier decomposition 
such as t he ones in^Z aritskv & Ri^ |l 99 7|): [B ourna ud et al 



van Evmeren et al 



(2005a|); lAngiras et all (|2006l . 12007^ 
(20119. 

To illustrate further. Figures Bl and B2 in Appendix 
B {electronic version) show some typical Hi maps from the 
WHISP sample for the minimum, mean and maximum val- 
ues of Asymmetry and Gini for all four lopsidedness cate- 
gories; none, weak, moderate and strong lopsidedness. In our 
view, these images illustrate how, for instance. Asymmetry 
and qualitative lopsidedness do not measure the same thing. 
A galaxy can be strongly asymmetric with a lot of fiux in 
a spiral arm offset from the centre of the galaxy, while at 
the same time, the outer contour may appear much like a 
ordinary disk. Conversely a strung-out galaxy may appear 
very lopsided at the lowest H I flux levels, but if there is 
little flux in the outermost part, and a strong, symmetric 
disk (with a ring for instance), this may not show in any 
of the morphological parameters. Our parameters are flux- 
weighted by design but the qualitative visual classification 
of lopsidedness by the previous authors may not be. 
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Figure 1. The normalized distribution of parameters of the combined samples from lSwaters et alj l|2002l ) and lNoordermeer et al.l l|2005bl ) 
for morphologically lopsided (dashed) and non-lopsided (dotted) histograms. There are small shifts in the distributions of Concentration, 
Gini and to a lesser extent M20 and Gm but no clear separation between the lopsided and non-lopsided galaxies. 



To verify if the shape of the outer contour alone is a 
better indication of lopsidedness, we compared the morpho- 
logical parameters for the images with uniform weighting 
(pixel- values set to h = 1). The Gini parameter, in this 
case, is of no use as this image is perfectly equal (G=0). 
The Asymmetry, M20 and Gm parameters show less change 
in distribution between the lopsided and non-lopsided pop- 
ulations. 

Alternatively, in order to parameterize morphological 
lopsidedness in similar terms as the morphological parame- 
ters presented, one could redefine Asymmetry using a spe- 



cific axis (requiring the additional input of a Position An- 
gle). For example, [Ba ldwin ct al. (1980) chose the east-west 
axis and lRichter fc Sancisil ||1994D chose the systemic veloc- 
ity axis. We compared fluxes from either side of the centre of 
these galaxies for the x and y-axes of the maps {A.j: and Ay, 
Tables Al and A2 in Appendix A, electronic edition) and 
found no relat i on wi th the lopsidedness qua l ifier fro m either 
ISwaters et all (|2002l) or lNoordermeer et all (|2005bh (see Ta- 
ble ?? and ??). Lopsidedness cannot easily be quantified us- 
ing the above common morphological parameters or simple 
variations thereof. At best, the histograms in Figure[l]can be 
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Figure 2. The distribution of parameters of the combined samples from lSwaters et all l|2002l ') and iNoordermeer et al for the 

different morphologically lopsided classifications (no, mild, moderate, and severe). There seems to be no clear preference of a lopsidedness 
classification for any part of the parameter space. 



used to assign a probability of morphological lopsidedness. 
In light of the fact that the above morphological parameters 
were developed to discern spirals from ellipticals in the opti- 
cal, their insensitivity to the lopsidedness of spiral H I disks 
is a indication of the extent that they can be used to clas- 
sify morphological s ub-types. T herefore, a Fourier decom- 
position (similar to 



l2005al : lAngiras et al 



Zaritskv fc R ix 1997; Bourn aud et al.l 
2006I . l2007i V is still needed to classify 



and quantify the lopsidedness of galaxies in future large H I 
surveys, such as WALLABY on ASKAP or a northern sky 



H I survey with APERTIF on WSRT or the MHONGOOSeQ 
nearby galaxy survey with the MeerKAT radio telescope. 



5.2 Interaction 

Based on the literature, one expects there be some signal of 
interaction in Concentration, Asymmetry, M20 and Gini as 



1 MeerKAT HI Observations of Nearby Galactic Objects; Ob- 
serving Southern Emitters 
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Table 1. The number of galaxies selected as interacting in the two WHISP subsets; the number and fraction of the sample, and number 
of individual galaxy in agreement with the visual classification by either author. The first three criteria are for the CAS and Gini/M20 
used in the literature. We defined criteria 3, 4 and 5 for H I morphology. The last three are various combination of our H I criteria. The 
criteria from the literature (1, 2 and 3) are for optical morphology and overselect Hi disks compared to the visual classification. Of the 
Hi criteria, (5) and (6) work well, with the latter agreeing with the visual classification in the case of individual galaxies. 





Noordermeer 




Swaters 






(68 galaxies) 




(73 galaxies) 




Selection Criterion 


Nr. (fraction) 


Individual 


Nr. (fraction) 


Individual 






agreement 




agreement 


Visual classification 


27 (39%) 




5 (7%) 




(1) A > 0.38 


55 (80%) 


23 


53 (78%) 


4 


(2) G > -0.133 X A/20 + 0.384 


51 (75%) 


20 


43 (63%) 


4 


(3) G > -0.4 X A + 0.66 


61 (81%) 


23 


65 (96%) 


5 


(4) Gm > 0.6 


39 (57%) 


21 


13 (19%) 


4 


(5) A < -0.2 X M20 +0.25 


22 (32%) 


2 


36 (52%) 


1 


(6) C> -5 X M20 + 3 


23 (33%) 


11 


8 (11%) 


4 


(4) & (5) 


6 (8%) 


1 


5(7%) 





(5) & (6) 


(0%) 





1(1%) 





(4) & (6) 


23 (33%) 


11 


8(11%) 


4 



the galaxy is warped and tidal arms are formed: the interac- 
tion spreads flux from the exponential distribution, altering 
Concentration and Gini, adds bright knots of stars further 
from the centre in optical images, changing both M20 and 
Asymmetry, and interaction breaks the overall symmetry of 
the galaxy's image, also modifying Asymmetry. Since this 
parameter space was developed to classify galaxy morphol- 
ogy, one can expect changes in more than one of the param- 
eters simultaneously when the spiral disk is gravitationally 
disturbed. During the interaction, one can also reasonably 
expect the disk to return periodically to unperturbed mor- 
phological values. Hence, we seek a section of this parameter 
space where interacting disks spend some of the several Gyr 
that an merger takes. The fraction of galaxies in this param- 
eter space, together with and estimate of the typical time 
spent there gives a typical merger rate for an H I survey. 

Figure [3] shows the histograms of the morphological 
parameters (Concentration, Asymmetry, Smoothness, Gini, 
M20 and Gm), for int eracting and non-in t eractin g, isolated 
galaxies, according to Noord ermeer et al.' (|2005bh . The five 
galaxies marked bv lSwaters e t aL (2002) as interacting are 
included. Again the interaction classification was done by 
a single observer, introducing some risk of a personal bias. 
The remainder of the Swaters et. al. sample is treated as iso- 
lated galaxies but it may contain some (mildly) interacting 
dwarf galaxies. In fact, given the observational result that 
lower mass galaxies at higher reds hifts (z=0.2-1.2) show high 



fractions of interactions (~ 10% Bridge et al.l 2007 
Kartaltepe et al.l l2007l: iLotz et alT 2008a : Lin et al 



201c 



2008 



Conselice et al.ll2009l : Ijogee et al.' 20091: iLotz et aLlbOlOal ). 
the real fraction of interacting galaxies may be higher. We 
therefore treat the Noordermeer sample as the cleanest and 
consider the Swaters sample more for confirmation. 

Both the CAS space and the Gini/M2o parameters have 
been used to identify morphologically disturbed galaxies in 
the literature and in the previous papers in this series, we 
established that the UV or FIR and the H I perspective trace 
similar structure. Thus, we expect to see some signal of in- 



teraction in many of these parameters, notably Asymme- 
try, M20 and Gm- The interacting galaxies show high values 
of Asymmetry ( A > 0.6). This i s somewhat higher than 
the cut used by IConselic3 (|2006l ) for optical asymmetry; 
Aopticai > 0.38 but they required A > S as well. The in- 
teracting galaxies have higher values of M20', M20 > -1, 
which is not too different from the cut used by iLotz et al] 
(2004). Direct cuts are most commonly used in these pa- 
rameter spaces. Alternatively, one could define the eigen- 
vectors of the interacting population in the combined pa- 
rameter space. However, simply excluding the locus of non- 
interacting galaxies and including most of the merging ones 
is the best one can do since the interaction qualifier is a 
subjective and qualitative one, not a q uantitative one like 
the tidal di sturbance parame t er used inlKarachentsev et al.l 
(200 41) and iBournaud erai] ll2005bl ). ISwaters et al.l (I2OO2I ) 
and iNoordermeer et alj ( 2005bl ) do not discern between 
weakly and strongly interacting. In addition, the training 
is likely too small to define eigenvectors. Therefore we use 
hard cuts in parameter space and compare how well these 
retrieve the fraction and objects that the Swaters and No- 
ordermeer papers marked as interacting. 

Following the example of lScarlata et al.l (|2007l ). we plot 
each morphological parameter against the others in Figured 
for the combined sample and for the Swaters and Noorder- 
meer samples separately in Figure [S] Both samples appear 
to occupy the same parameter space, so combining them 
is not an issue. Interactin g galaxies are m arke d. Figure [5] 



confirms the assertion from IConseUcd l|2003t ) and lLotz et al.l 
(|2004h that Asymmetry and M20 are parameters sensitive 
to m ergers, and our assertion in .Holw crda et al. (2011b:3, 
that Hi is a good wavelength to investigate it. We 
define three criteria to select mergers: 



Gm > 0.6, 

A < -0.2 X M20 + 0.25, 
and 



(9) 
(10) 
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Figure 3. The distribution of parameters of the combined samples from ISwaters et al and iNoordermeer et al for 

interacting (dashed) and non-interacting (dotted histograms). There are clear separations in the distribution of Asymmetry, M20 and 
Gm values between the two populations of galaxies. Values in Tables Al and A2 in the electronic version of the manuscript. 



C> -5 X M20 + 3. 



(11) 



G > -0.4 X A + 0.66 or A> 0.4 



(14) 



These are the dotted lines in Figures |4] and [5] (equation [9] in 
Panels I, III, VI and X, equation 1101 in panel V and equa- 
tion [11] in panel IX) . We also use three criteria for these 
parameters defined in the Uterature: 



A > 0.38 

G > -0.115 X M20 +0.384 
and 



(12) 
(13) 



With the first one fro m IConselicd l|2003l) and the last two 
from iLotz et all (|2004 l2010bl ). These are the dashed lines 
in Figures |4] and [5] panels IV- VI, II and IV respectively. 

Combined with those from the literature, we list their 
success rates in Table [1] From this table, it is evident that 
the criteria from the literature do not translate well to H I 
column density map morphology. These criteria select too 
many contaminants. In part this may be because both Gini 
and Concentration are linked and sensitive to how concen- 
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Morphology parameterspace 
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Figure 4. The parameter space of the combined samples from ISwaters et al and iNoordermeer et al ] (2005b') for interacting 

(black symbols) and non-interacting (gray symbols). There are clear separations in Asymmetry, M20 and Gm distributions between the 
two populations of galaxies. Our two cuts in parameters are indicated with a dotted line; A > 4.3M20 and Gm > 0.6. Combined, these 
two cuts select a reasonable fraction of the interacting galaxies. Values in Tables Al and A2 in the electronic version of the manuscript. 



trated the image is. As H I maps a re more extended, ther e 
is a shift in values (see paper II, iHolwerda et~all l201ld ). 
The Gm parameter criterion performs well, selecting most 
interacting galaxies with Gm > 0.6 but with quite some 
contamination, even in the Noordermeer sample (Table [T|) 
which is the cleanest of the two. Therefore, a combination 
of one or more morphological parameters appears the most 
promising to cleanly separate an H I sample into interacting 
and isolated galaxies. The Asymetry-M2o selection criterion 
performs better in that it selects a similar fraction of galax- 



ies but it does not agree with either Swaters or Noorder- 
meer estimate in the case of individual objects. It does se- 
lect many more objects in the Swaters sample but as we have 
noted above, one can reasonable expect mor e dwarfs to be 
(mildly) interacting than the five flagged by ISwaters et al.l 
(2002). The Concentration- M20 criterion works best as it se- 
lects a similar fraction of galaxies as the visual classification 
but also agrees on more cases of individual galaxies. It also 
agrees well with the Swaters selection of interactions. Al- 
ternatively, a Concentration- Asymmetry criterion may well 
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work. For instance, combining any two of these criteria does 
not improve the selection appreciably. Combining the Gm 
criterion and the C/M20 criterion, effectively is the latter 
criterion. 

We intend to apply these morphological cuts on rep- 
resentative samples of Hi observations, starting with the 
complete WHISP sample (|lfolwerda et aLllioilal ). To con- 
vert these fraction into a volume merger rate, one needs 
to compute the representative volume of the survey and a 
timescale for which merging systems reside in the interaction 
part of parameter space. We focus on these time scales in the 
next paper in this series (|Holwerda et al.ll2011dh , using SPH 
simulations of gas-rich 1:1 mergers. 



6 CONCLUSIONS 

Based on ou r quantified morphologic al analysis of 141 
jalaxie s from ISwaters et al., (2Q02) and iNoordermeer et al.l 
2005bl ) for which they provided visual estimates of the lop- 
sidedness and level of interaction for Hi disks, we can con- 
clude the following: 

1. The two-dimensional morphological parameters can- 
not discriminate between weak and strong lopsidedness as 
judged visually by previous authors. However, Asymmetry, 
M20 and Gm , and to a lesser extent Gini parameters all show 
a shift in the mean of the distribution of values between the 
lopsided galaxies and those that are not lopsided (Figure O. 

2. We suggest, therefore, that these parameters can be 
used to assign a probability of lopsidedness (Figure [T]). But 
future surveys should use Fourier analysis to find lopsided- 
ness in the Hi distribution. 

3. The fraction of interactions in a sample of H I maps can 
however be determined similarly well using these parame- 
ters, as a visual classification. Individual parameters, such 
as Asymmetry and Gm, do not select the interacting sys- 
tems cleanly (Figure [3] and |4|. 

4. Combined criteria, using Asymmetry and AI20 or Con- 
centration and M20, work better (Figure U and [SI and Table 
[1]) and select the right fraction of a sample of galaxies is 
currently undergoing interaction, as identified by visual in- 
spection. Combined with an estimate of the time a merger is 
selected by these criteria, one can estimate what the merger 
rate in an Hi survey is. The benefits of such a merger rate 
determination would be less observer bias than a visual clas- 
sification and a empirical visibility time, determined from 
simulations. 



7 FUTURE APPLICATIONS 

The parameter space, as we applied it to the WHISP H I 
column density maps, allows us to find candidates for lop- 
sidedness and more accurately define the fraction of interact- 
ing galaxies, solely from their H I morphology. It remains to 
be determined how long an interacting disk remains in the 
interaction part of the morphology parameter space. This 
can be addressed with the new generation of simulations of 
major and minor mergers currently being u ndertaken (e.g. 



Bournaud et al 



r merg e rs currently being u ndertaken (e.g, 
l2005bl: ICox et all [2006al lbl: IWeniger et al] 



20091 : iLotz et al.|[2010al lbl). which include a comprehensive 
treatment of the interstellar matter in the galaxies during 



the merger. The time-scale for which a disk has a morpho- 
logical interaction signature can then be determined by aver- 
aging over the many possible viewing angles. This timescale 
and the full WHISP sample (368 galaxies) will allow us to 
estimate the interaction rate of spirals locally, based purely 
on their Hi morphology. This can serve as an additional 
zero-point for estimates of the merger rate at higher red- 
shift. Upcoming nearby galaxy surveys with MeerKAT and 
the Widefield ASKAP L-band Legacy All-sky Blind surveY 
(WALLABY, Koribalski et al. in preparation) for the South- 
ern Sky and Northern Sky Survey with APERTIF on WSRT 
will solidify the local Universe merger rate estimate, based 
on H I morphology. The future Square Kilometre Array can 
subsequently determine the merger rate of gas-rich galaxies 
over Cosmic times (up to z~l or better). The great benefit 
of H I surveys to determine the merger rates are the sensitiv- 
ity of H I to interaction and the sensitivity of H I surveys to 
lower mass syste ms, for which the merger rate is the poorest 
constrained (see lLotz et al.ll2010al lbh. 
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